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Abstract 

We study the long-time behaviour and the spatial correlations of a simple ferromagnetic spin 
system whose kinetics is governed by a thermal bath with a time-dependent temperature which is 
characterized by a given external relaxation time r. Exact results are obtained in the framework 
of the spherical model in d dimensions. In the paramagnetic phase, the long-time kinetics is shown 
to depend crucially on the ratio between r and the internal equilibration time Tcq. 

If r < Tgq, the model relaxes rapidly towards an equilibrium state but there appear transient 
and spatially oscillating contributions in the spin-spin correlation function. On the other hand, if 
r S> Teq the system is clamped and its time evolution may be delayed with respect to the one of 
the heat bath. For waiting times s such that r ^ s ^ Teq, a quasi-stationary state is found where 
the fluctuation-dissipation theorem does not hold. 



^PACS numbers: 05.70.Ln, 05.70. Jk, 05.40+j, 82.20.M 



^Laboratoire associe au CNRS UMR 7556 
^Laboratoire associe au CNRS UMR 7085 



1 Introduction 



Non-equilibrium systems generally display specific features which are absent from their corresponding 
thermodynamically equilibrated counterparts. For example, consider a many-body system contained in 
a thermostat and prepared in some initial state. Then assume that the temperature of the heat bath is 
rapidly quenched to a different temperature Ti and the system is then allowed to evolve freely, with the 
temperature fixed at Ti. In this setting, the system under study may show ageing effects, which come 
from the dynamical breaking of time-translation invariance, see 13 El El for reviews. Although ageing 
is known, since a very long time, to occur in glassy systems, similar effects have more recently also been 
observed for non-disordered spin systems. In simple ferromagnets such as the Ising model, ageing effects 
have been studied through the long-time behaviour of two-time correlation and response functions, see 
jH] and references therein. In addition, it has been proposed recently that the well-known dynamical 
scale invariance, associated with the ageing phenomenon (see e.g. 5112]) in simple ferromagnets, may 
be generalized towards a larger group of local dynamical scale transformations ^ . Among other results 
this theory leads to an explicit prediction for the time-dependent thermoremanent magnetization which 
has been confirmed in several distinct models [7| 1511^ ITUIITT] . 

However, current theoretical studies generally assume that the quench is infinitely rapid, or in other 
words that the final temperature Ti is reached immediately. In real experiments, this condition is often 
far from being satisfied, see e.g. jT^i El- It is therefore interesting to investigate the behaviour of a 
system coupled to a bath with a time-dependent temperature T{t). Here, we shall study the long-time 
behaviour of a system initially prepared in a fully disordered state and brought at time t = into 
contact with a thermal bath which evolves from an initial temperature Tq to a final temperature Ti 
in a finite or infinite interval of time. In the present paper, we shall study the following form of a 
time-dependent temperature 



where r is a given relaxation time which sets an external time scale Text = t- Certainly, the long-time 
behaviour of a statistical system will depend, besides on the external time scale r, also on the internal 
time scale(s) rint(T) which arise from the internal self-organization of the system, governed by the 
microscopic interactions. Two basic cases should be expected. First, if r^^t ^ t, the system simply 
relaxes towards equilibrium at the final temperature Ti. On the hand, if Tint <^ t, the system is said 
to be damped at the external temperature T{t) and its time history will depend on the properties of 
the external heat bath. While for times t ^ r ^ Tint the system should have simply relaxed back to 
equilibrium, new effects might be expected to occur in the regime r S> t S> Tint and this may allow to 
distinguish between these two cases. The time-dependence as specified in eq. flLlj) should contain the 
basic ingredients which are needed for a conceptual understanding of this kind of phenomenon.^ 

One of the central questions is whether/when under the conditions just described the system may 
evolve towards thermodynamic equilibrium. A convenient way to measure the distance of a system 
from equilibrium is through the fluctuation-dissipation ratio [1^1 ITB] 



where C{t,s) and R{t,s) are the two-time auto-correlation and autoresponse functions, respectively 
(see section 2 for the precise definitions). At equilibrium, the fluctuation-dissipation theorem states 

■^Originally, we became interested in this problem through the consideration of the fragmentation of microscopic systems 
like atomic nuclei where it has often been assumed that the detected fragments are at thermodynamic equilibrium and 
fixed temperature ^3] . However, such descriptions do not take care of the time evolution of these systems which expand 
in space over a finite time interval and cool down before they are detected. 
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that X{t,s) = 1. Therefore, it is an important question under which physical conditions equihbrium 
can be reached. This question has received much attention in glasses, see |17| ITH] for recent theoretical 
and [T^ IT^ 120] for recent experimental examples. Understanding this point must come before the more 
complicated questions relating to an eventual ageing behaviour can be addressed. 

In order to obtain explicit analytical results, we shall study the exactly solvable kinetic spherical 
model with a time-dependent bath temperature, to be defined precisely in section 2. For constant 
temperatures, this model has been studied in great detail in the past, either in the context of continuum 
field theories EH 1221 123 121123 13 or else in the form of a lattice model [2nil23l2Hllinil2nilHlEn]- It 
is well-established that the spherical model, despite its technical simplicity, still contains the main 
physical features encountered in other systems (such as the Ising model) which might be considered 
to be closer to the experimental reality but which no longer permit an exact analytical solution. In 
particular, results found from the spherical model in d < 4 space dimensions are known to be different 
from the predictions of mean-field theory. In this paper, we shall be concerned exclusively with the 
question under what condition a stationary or quasi-stationary state of the model may be considered 
to be an equilibrium state in the sense that X(t,s) = 1. This question can be studied by restricting 
attention to the high-temperature phase such that the system is brought into contact with a heat bath 
at the initial temperature Tq. The heat bath then cools down to the final temperature Ti such that 
both are above the critical temperature, viz. Tq > Ti > > and the model remains in the disordered 
high-temperature phase. In that phase, there is no ageing behaviour to be expected. We leave the 
question of a possible ageing behaviour in systems coupled to a time-dependent bath for future work. 

The content of the paper is as follows. In section 2 we introduce the model and obtain the formal 
exact solution. The equal-time correlator is analysed in section 3 and we discuss several new non-trivial 
time and length scales which arise. Two-time quantities and the fiuctuation-dissipation ratio are studied 
in section 4 and the results are generalized to include the clamped case in section 5. Section 6 presents 
our conclusions. 



2 Model and formalism 

We begin by recalling the definition of the kinetic spherical model, using the formalism as exposed in 
|2nil2HllH]- We consider a system of time-dependent classical spin variables Sx(t) located on the sites x 
of a (i-dimensional hypercubic lattice. They may take arbitrary real values subject only to the spherical 
constraint 

J2s.itf=^ (2-1) 

X 

where A/" is the number of sites of the lattice. The spherical model Hamiltonian reads 

n=-J J2 Sx{t)Sy{t) (2.2) 

<x,y> 

where the sum extends over nearest-neighbour pairs only. In the following we choose units such that 
J = 1. The system is supposed to be translation- invariant in all directions. The kinetics is assumed to 
be described in terms of a Langevin equation 



dt ^ 



Sy{t)-{2d + i{t))Sx{t)+r^x{t) (2.3) 



where the sum over y extends over the nearest neighbours of x and rix{t) corresponds to a stochastic 
force which describes the action of an environment which lies outside of the nearest-neighbour range. 
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Physically, this means that the model is assumed to be immersed into a heat bath. The resulting forces 
are supposed to be gaussian and thus to be characterized by an ensemble average which is zero, viz. 
(Vxit)) = and a second moment which reads 



{Vxit)Vy{t')) = 2T{t) S^Jit - t') 



(2.4) 



At time t = the system is brought into contact with the heat bath. We assume in this work that the 
temperature of the heat bath decreases from To to Ti with a characteristic relaxation time r according 
to 

T{t) = Ti + (To - Ti) exp(-t/r) (2.5) 

The choice of this temporal behaviour is dictated by computational simplicity. In addition, as we shall 
see, it already contains the basic effects which will also be present for more general thermal histories. 
For sufficiently large values of the initial bath temperature Tq, the typical correlation lengths are of the 
order of a lattice constant or less. Then for all practical purposes, the system is effectively uncorrelated. 
Finally, the function lit) is fixed by the spherical constraint ()2.H) and has to be determined. In order 
to do this, we follow the standard procedure of replacing the spherical constraint by its mean value. It 
can be shown that this does not change results if the infinite-system limit A/" — >■ oo is taken before the 
long-time limit t — oo [30j. 

By a Fourier transformation 



(2.6) 



where the integral is taken over the first Brillouin zone E, the Fourier-transformed spin variable S{q, t) 
becomes 



S{q,t) 

with the dispersion relation 
and we have also defined 



-uj(q)t 



UJ 



S{q,0)+ / dt'e-('^)*'v^^(q,t') 



;q) = 2 5^(l-cos(g.)) 



i=l 



g{t)=exp(^2 j^dt'iit'] 



(2.7) 



(2.^ 



(2.9) 



Clearly, the time- dependence of S{q,t) and any correlators will be given in terms of the function 
9 = 9{t). 

We now derive the expressions for the correlators and response functions for an arbitrary thermal 
history T = T{t). We begin with the equal-time spin-spin correlation function 



C^,y{t) = C^_y{t) = {S^{t)Sy{t)) 

We obtain the Fourier transform C{q,t) from 

{S{q,t)S{q',t')) = {2n)''6{q + q') C{q,t) 

and immediately find 



Ciq,t) 



git) 



C{q,0) + 2 / dt' T{t') e''^^''^' g{t' 



(2.10) 



(2.11) 



(2.12) 
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If we assume uncorrelated initial conditions^ 

Cx,y{0) = 5cc-y,0 (2-13) 

as we shall always do in the following, we have C{q,0) = 1. Because of the spherical constraint (|2.2j) 
and spatial translation invariance, the autocorrelator must satisfy 

Coit) = [ dqCiq,t) = (SUt?) = 1 (2.14) 
Jb 

This in turn fixes 3(t) or via (|2.9|) the function g{t) as the solution of a Volterra integral equation 

g{t) = fit) + 2 fdt' Tit') fit - t')git') (2.15) 
Jo 

where 

fit) = ^ jdq e-2-('^)* = {e~''hiAt)Y (2.16) 

and Jo is a modified Bessel function jSI]. Once we have found the solution of eq. ()2.15|) . the correlation 
function can be obtained. 

Before doing this, we now give the expression for the two-time correlation function Cx-yit,s) = 
{Sxit)Syis)) and the two-time response function Rxit,s). The calculation follows entirely standard 
lines ||26||2H1E] and we merely quote the result. In Fourier space 



C(g,t,.)=C(g,.)e-('')(*-^)y^ 



(2.17) 



and where g = git) is the solution of eq. ()2.15p . Similarly, the response function is obtained in the usual 
way [231211123123111 by adding a small magnetic field term 6H = — J2x hxit)Sxit) to the Hamiltonian. 
We easily find in Fourier space 



5hiq,s) 



hr=0 



^-.i,)it-s) 9i^ (2.18) 

^' git) 



From these expressions, the autocorrelation function C(t, s) = Coit,s) and the autoresponse function 
i?(t, s) = Roit, s) can be obtained by integrating over the momentum q. 

Summarising, the physically interesting correlation and response functions are given by equations 
()2. 1212. 1712. 18|) together with the constraint eq. ()2.15|) . for any time-dependent temperature T = T(t). 
This constitutes the main result of the general formalism. 

In order to solve eq. ()2.15p explicitly, we now use the specific form eq. ()2.5p of the time-dependent 
temperature T = T(t). Through a Laplace transformation 

POO 

Jip)= / dtfit)e^^' (2.19) 
Jo 

eq. ()2.15p is transformed into a linear difference equation 

9ip)=7ip) + 2T Jip)gip) + 2(To - T,)Jip)gip + 1/r) (2.20) 

^The treatment of correlated initial conditions, following the lines of j8j, presents no additional difficulty. 
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The analytic solution of this equation allows to study the behaviour of the correlation and response 
functions defined above. Two limit cases correspond to constant temperatures and have been analysed 
in detail in the literature [2H]- First, in the limit r oo, we have a constant temperature T = Tq and 
the final temperature Ti is never reached. Second, the limit r — corresponds to an infinitely rapid 
quench to the constant end temperature T = Ti. We are interested in the behaviour between these two 
extremes. The solution 'g{p) of eq. (|2.2(Jj) can be cast in the form 



m = Y.^T,-T,rup) (2.21) 



n=0 

where 



^n(p) = 2"n^(°)(Ti, p+^'j (2.22) 
fe=o ^ ^ ^ 

Here 'g'^'^^ {Ti,p) is the solution as given by Godreche and Luck for a fixed temperature T = Ti 

^(o)(T,,p) = (2.23) 

We shall turn in the next section to a detailed analysis of the properties of the physical observables 
coming from this solution. 



3 The equal-time correlator 

We now analyse the long-time behaviour of the equal-time spin-spin correlation function Cx{t) or rather 
its Fourier transform C{q,t). Throughout this paper, we shall restrict ourselves to the situation where 
both the initial temperature Tq and the final temperature Ti are above the critical point and the systems 
is cooled from Tq to Ti, viz. 

To > Ti > T, (3.1) 

The equilibrium critical temperature is given by Tc = (2/(0))~^, see [23123123; and is non-zero for 
d > 2. As we have seen in the last section, the time-dependence of the correlators follows from the 
form of the function g{t). In turn, the behaviour of g{t) can be described in terms of the singularities 
of its Laplace transform 'g{p). From eqs. ()2.21I2.22|) . these singularities are entirely given in terms of 
the singularities of the constant-temperature solution 'g^^\Ti,p) of eq. ()2.23|) . 

In the high-temperature phase, the singularities of g^^^ (Ti , t) can be analysed as follows j2H] ■ Consider 
the function 

which is monotonously decreasing with p. Furthermore, f{p) is analytic in the complex p-plane except 
for a cut in the interval —8d < p < 0. On the other hand, if Ti > T^, the function g^^''{Ti,p) has 
a simple pole at a value po given by /(po) = l/(2Ti). This translates into an exponential long-time 
behaviour of g^^\Ti,t) ~ exp(t/req) which defines an equilibration time 

Teq := req(Ti) = Tq = — (3.3) 

Po 
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Now, for a time-dependent temperature T = T{t), the only possible singularities of g{p) are those 
of 'g^^\Ti,p + k/r) and therefore the singularities of g{p) with a positive real part are simple poles and 
occur at 

k 

Pk ■=Po (3.4) 

T 

with A; = 0, 1, 2, ... , fcmax. Here fcmax is the largest integer such that Pk^^^ is still positive. For k > k^ax, 
the singularities merely lead to exponentially decreasing corrections in g{t) and do not contribute to 
the leading long-time behaviour we are interested in. The sequence of poles pk leads to the following 
sequence of relaxation times 

1 _ 1 k 

The quantities pk are positive for < A; < k^^^y^ and negative for k > fcmax. Hence for t sufficiently large, 
i.e. t > req(Ti) the leading contributions to g{t) are of the form 

^max 

g{t) ^ 5^ Gfe exp(t/rfc) (3.6) 



k=0 



where 



and in addition 



j=0,j^k 



T 



n=0 



7 = -/(Po)7/ (Po) 
is a positive constant and the prime denotes the derivative. 



(3.7) 
(3.8) 



Using the expression eq. ()3.6|) the spin-spin correlation function given by eq. ()2.12p takes the following 
form 



C{q,t) c^T^^J2Gkexp{-kt/T) 



,fc=0 



Gn 



-ft/ 



Gk 



1 '^fcmax'f'l 



,-(A:max + l)tA 



uj{q) + A-q2 uj{q) + 

provided t > Teq(Ti) and where for all < £ < fcmax + 1 we have defined 



(3.9) 



2req(Ti)r 



r - £req(Ti, 



1/2 



(3.10) 



In the following, we shall use the relation Aeq = Aq = \/^2rr^. In addition, we also used the abbreviation 

(3.11) 



Hk — Gk -\ ^^t;; Gk-i 



In order to understand the behaviour of G{q,t), we note that the first term in eq. ()3.9|) . which 
formally corresponds to £ = 0, gives the stationary contribution. In addition, there is at least one 
additional transient term which will disappear in the limit t ^ t. The number of transient terms which 
are present depends on the value of fcmax = fcmax(^i5 t)- We now discuss the various physical cases which 
can arise. 
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3.1 First-mode contribution — Short thermal relaxation times 



The simplest case arises if T^q > r, that is the bath relaxes faster than the system itself. Thenpo~ ^ 0, 
that is /cmax = and we have in eq. ()3.9p £ = as the only contributing mode. Therefore (assuming 

t > T and t ^ Tcq(Ti)) 

ciq, t) ^ T, ( ,^\_,,^, + + ^ i^-''^^'^'-"') (3.12) 

For a physical understanding, it might be more appealing to consider the correlator in real space. In 
the ID case, the result is particularly simple and becomes in the large-separation limit |a;| ^ 1 

C^{t) ~ Ti (^Aoe"l^l/^« - ^^^Ai sin (^M^ g^*/^^ + O (e-^'/^e"*/^) (3.13) 

Therefore, the correlator is the sum of a stationary and spatially homogeneous term and a transient 
term which in addition shows spatial oscillations of wavelength Ai. We point out that Ai diverges when 
the external time scale r approaches the internal equilibration time Teq from below, see eq. (j3.1(jp . In 
figure H we show the location of the line Teq = r as the curve 1/ri = in the (Ti, r)-plane. The region 
I in figure HI corresponds to the case fcmax = 0. 

For larger values of r we go over to a multimode regime which we discuss below. 




r I 

1 



T 

Figure 1: Kinetic phase diagram of the 3D spherical model with a time-dependent temperature as 
given by eq. (j2.5j) . The lines shown are the loci where the relaxation times r^, £ = 1, 2, 3 diverge. The 
dotted line marked gives the equilibrium critical point and the regions I, II and III are those with 
kuiax = 0, 1, 2, respectively. 
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3.2 Multimode contributions 



If ^max > several transient modes contribute to C{q,t). Their contributions are quite different, 
however. The propagator arising in all terms is of the form 

1 q^O 1 



uj{q) ± 



(3.14) 



with a positive sign for k < k^ax and a negative sign for k = fcmax + 1- This leads to a different spatial 
behaviour, in analogy with eq. (j3.13|) found for the case /Cmax = discussed above. The first k^ax terms 
all lead to spatially decaying contributions to Cx{t) in the large- 1 a; | limit, where the characteristic length 
scale is given by Xe with £ = 0,1,..., fcmax- With the exception of the i = contribution, all terms 
are transient and decay with a relaxation time given by r/£. In contrast, the last term, which is also 
transient, gives rises to a spatially oscillating term with a wave length Aosc := Xk^^^+i- 

The behaviour of the system in the first few regimes is illustrated in figure H For a fixed final 
temperature Ti and increasing values of the external relaxation time r, the system goes over from 
region I with only a single mode (since k^^ax = 0) to region II with modes, later to region III with three 
modes and so on. 









I ' ' ) 




\ / 




\ 1 




\ / 




\x, / 

\ 1 / 

\ / 

\ / 

\ / 
\ / 




vy 

/ \ 
/ \ 
/ \ 




/ \ 






/ 








/ I 


II 

1 


i 


III 

1 







0.5 



Figure 2: Dependence of some correlation and wave lengths on the external relaxation time r, at fixed 
final temperature Ti = 3.8363 ... in the 3D spherical model. The regions I, II, III correspond to the 
cases fcjnax = 0, 1,2, respectively. The thick solid lines give the oscillation wave length Aqsc = -^fcmax+i- 
The dash-dotted lines give the equilibration length scale Acq = Aq (the equilibration time Teq ~ Agq) 
and the first additional length scales Ai and A2. The dashed hne corresponds to the curve A = \/2t, see 
text. 

The behaviour of the various correlation and wave lengths is illustrated in figure |21 for the 31^ case 
and for T^q = 0.4, that is Ti = 3.8363.... First, we note that the characteristic length scale Aeq 
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associated with the only non-transient term does not depend on the externally imposed time scale r. 
Second, the thick solid lines show the dependence of the oscillation wave length Aosc = -^fcmax+i If 
we start with a small value of r in region I and then increase r, Aosc increases and finally diverges when 
T = Teq that is whcn the external time scale becomes equal to the internal equilibration time scale. If r 
is now increased further, we enter into region II. At this point, the previous oscillating mode transforms 
into a spatially monotonous one which simply gives another contribution with decreases as a function 
on I a; I and in addition, a new oscillating mode with a finite wave length appears. This new mode arise 
since one of the modes which in region I leads to exponentially decaying contribution to g{t) transforms 
itself into an, albeit non-leading, exponentially growing contribution. Furthermore, the characteristic 
length of the old oscillating mode (denoted by Ai in figure 12]) now decreases as a function of r. We now 
increase r further until the oscillation wave length Aosc diverges again and we therefore go over from 
region II into region III. At this transition point, we observe that the value of Ai agrees with the one 
of the new oscillation wave length for region III. Phenomenologically, this looks as if that mode would 
split into two modes, one giving rise to a spatially oscillating contribution to Cx{t) and the other one 
to a spatially decaying term. This splitting occurs along the line A = \/2t and in particular, it follows 
that the oscillation length always satisfies Aosc > v^2r. 



4 Two-time correlations and response functions 

Using the expression of the function g{t) from eq. (j3.6|) . we can formally compute the two-time ob- 
servables we are interested in, that is the two-time autocorrelation function and the autoresponse 
function defined in eqs. ()2.17l2.18p . In particular, a study of these two functions allows to calculate 
the fluctuation-dissipation ratio X{t, s). A deviation of that ratio from unity should signal a departure 
from thermodynamic equilibrium. For a sudden quench to a final temperature Ti this question has 
already been studied j2H] and it was shown that indeed X{t, s) = 1 in the long-time limit in the high- 
temperature phase when Ti > T^. Here, with the additional complication of a time- dependent heat 
bath temperature, we must first define an instantaneous time-dependent temperature of the system in 
order to make the criterion of the fluctuation-dissipation ratio applicable. We choose the equilibrium 
part Ceq of the single-time correlator as the physical quantity which measures a temperature. Consider 
for a moment the case of a fixed temperature T. The equilibrium part behaves for small momenta (or 
large distances) as 

limaq(9) = lim — = 2reqT (4.1) 

and depends on T in a known way. Therefore, we use the extension of this relation to define, for a 
time-dependent bath temperature T{t) varying according to eq. ()2.5|) . an effective temperature Tefj(s) 
of the spin system. Taking the explicit form ()3.9|) of the single-time propagator C(g, s) we set 

limC(q,s) = 2r,q(Ti)Teff(s) (4.2) 

q^O 

which is the defining equation for the effective time-dependent temperature Tcff(s). It remains to be 
seen whether this Teg can play the role of a physical equilibrium temperature. 

This can be seen by means of a comparison with the results obtained from an instantaneous quench 
towards Ti, hence r = 0. For clarity, we shall denote by A^^'^ the result obtained for an observable A, after 
a sudden quench with r = to the final temperature Ti. The two-time autocorrelation function C^'^\t, s) 
is explicitly given in eq. (2.58) of |2H!. In particular, it is known that [2HI liniq^o C^^Kq, s) = 2req(Ti)ri. 
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We now consider the two-time autocorrelation function C{t, s). From section 2, we have 



9{t) 



44 /dg^CW(0,.)e-'^^(*-^) 
9{t) J Ti 

k^^\t,s)C^''\t,s) (4.3) 



Ti 

In the second and the fifth fine we used the fact that we are in the regime t — s ^ 1 , in combination 
with the fact that in this hmit, only the low-q behaviour of C{t, s) and of C^^\t, s) really contributes. 
In the fourth line, the definition of Tes{s) was used. For brevity, we have set 



where the Laplace transform of g^'^\t) was given in ()2.23|) . The exact autoresponse function is given by 



R{t, s) = Jdq Riq, t, s) = f (^^^ \J = ^^"^^ 

We now ask whether the quasi-equilibrium temperature Tefj(s) so defined is really the temperature 
of a system in equilibrium. A convenient way to test this is through the fiuctuation-dissipation ratio 

X{t, s) = ns{s)R{t, s) (^%^) (4.6) 

An equilibrium situation can only be recovered if X{t,s) = 1. We now bring this expression into a 
simpler form as follows 

1 dC{t,s)/ds 



X{t,s) T,s{s)R{t,s) 

1 Cit,s)dCit,s)/ds 
%s{s)R{t,s) C{t,s) 

1 CW(t,g)reff(g)Tf^fcW(t,g) d 

reff(s) R(^\t,s)m{t,s) ds 

1 C(o)(t,s) d 



lnC(t,s) 



TiRio){t,s) ds 



\nCit,s) 



9j^C{t,s)j (^-lnC(°)(t,.)) (4.7) 
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In performing this calculation, we used in the third line the approximate expressions eq. (j4 .314.51) valid 
for t — s ^ Teq{Ti) and in the last line, we used the fact the fluctuation-dissipation theorem holds for 
any instantaneous quench to the final temperature Ti > Tc j2H]- 

After some algebra, we can write this in a compact form, taking also the t — >■ oo limit 



X{t,s) 



1 + 



T 



(4.8) 



Remarkably, this scaling form only depends on the single scaling variable x = s/r > 1. The associated 
scaling function reads 

d d 
^(x) = ( — In $i(a;) + 2— In $2(3^) 



where explicitly 



dx 



dx 



(4.9) 



$i(a;) 



$2(0;) 



(-TO 



-kx 



k=0 



Ti 



.e=o 



-ix 



Go^lq H 7^ — Gk^^Xl^ 

J 1 



^ + J2 HkXle- 



kx 



k=l 



(4.10) 



(4.11) 



A particularly simple form is found for A;max = 0. Physically, this corresponds to the situation when 
T < Teq(Ti). Then the fluctuation dissipation ratio reads 

X{t,s) = 

As in the case of an infinitely rapid quench, the fluctuation-dissipation ratio is unity, up to corrections 
which vanish exponentially fast in the long-time limit. We can conclude that the system evolves towards 
thermodynamic equilibrium. 



'eg 
r 



A?(To-Ti: 



■ exp(x) — 1 



X 



s 

T 



(4.12) 



5 The clamped case and the limit T/Teq(To) oo 

Studying the kinetics with a time-dependent temperature and in the disordered phase, way may identify 
the following time regimes: 

1. r ^ roq(Ti). This case formally corresponds to k^ax = 0. Studying the long-time kinetics, we 
always work in the situation where t 3> req(Ti) ^ r. Up to exponentially small corrections, we 
recover an equilibrium dynamics as found for an infinitely rapid quench [2H]- 

2. r ~ rcq(Ti). In this case, we consider t ^ req(To) and also t ^ r. This regime corresponds to a 
finite non-zero value of fcmax- With respect to the situation of an infinitely rapid quench, there are 
several new exponentially small contributions which arise, as was shown in section 4. As expected, 
equilibrium dynamics is recovered. 
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3. r ^ req(Ti). This is called the clamped regime^ because the microscopic degrees of freedom 
follow the varying temperature of the external heat bath. Certainly, for times t ^ r ^ Teq(To) 
we shall simply recover equilibrium dynamics, along the lines of section 4. Because we con- 
sider only quenches which reduce the temperature Tq > Ti > throughout, we always have 
Teq(^i) > TeqiTo). On the othcr hand, new and interesting long-time behaviour occurs in the 
regime Teq(To) ^ t ^ r. In this case, the system is in equilibrium as far as the internal degrees 
of freedom are concerned, since t ^ req(To), but has not yet equilibrated with respect to the 
external time scale of the heat bath, since t <^ r. (a special case of the case under consideration 
here corresponds to /cmax = 00.) 

We first have to look closer at the limit r — > 00. As we have seen above, the central part of the 
calculation is the solution of the Volterra integral equation ()2.15j) coming from the spherical constraint. 
However, for r — > cxd and after having performed the Laplace transformation, we may expand in 1/r 
and then obtain to leading order an ordinary differential equation 

^ + raip)g{p) = rb{p) (5.1) 

where _ 

a{p) = '^'^^P^-J , b{p) = -i— ^ (5.2) 
^"^^ 2(ro-Ti)/(p) ' ^"^^ 2T0-T1 ^ ^ 

Although this equation may be solved straightforwardly through integrations, it is for our purposes 
more useful to proceed differently. After all, we are merely interested in the large-time (or small-p) 
behaviour of the correlators. In addition, because of the r — > cxd limit implicit in eq. dOJ, we must 
take the limit in such a way that Ta{p) remains finite. Therefore, we merely need a{p) in the vicinity 
of pc defined by a{pc) = 0.^ The leading long-time behaviour should then be given by an expansion in 
powers of p — Pc- We write, using ()5.2|) 



. X A . N A df{p) 

a{p) -A{p-p,) + ... , A- ° '^"^^ 



To - Ti dp 



< (5.3) 

P=Pc 



We look for a solution of 1)5.11) of the form 

^(P) = E^~"^"(P) (5.4) 



n=0 



and we directly obtain 



In order to perform the inverse Laplace transformation, we write 



n=0 ^ n=0 ^ ^ 



e-(P-P^)* (5.6) 



^en frangais: serre; auf deutsch: cingcklcmmt 

^Physically, one may understand this by noting that in the case of constant temperature Tq, t is. formally infinite and 
we should have to leading order g{To,p) — b{p)/a{p). Therefore, the root pc sets the entire physical behaviour, up to 
normalization factors. 
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At this stage the sum may be performed first and we find 



oo 




9{t) = ^gn{t)e 



(5.7) 



ra=0 



which is vahd provided t <^ t. In this expression, Tefj is the typical time after which the function g 
and then the correlators which are closely related to it differ notably from the expressions they take at 
formally infinite external relaxation time r. We find 



Remarkably enough, Tefj differs from the external time Text = t- Rather, it is the geometric mean between 
the externally imposed relaxation time scale as measured by r and an internal time scale, measured 
by Tdep = 2A(To,Ti). In the physical observables, t^q appears as a natural time scale informing on the 
departure from an effectively constant temperature T = Tq. We point out that this new time scale Tdep 
is distinct from the equilibration time at temperature Tq, namely req(To) = l/pc- 

For example, the two-time response function becomes 



which illustrates again the role of t^s in the description of the deviation with respect to R^^\t, s). 

We now comment on how these results depend on spatial dimensionality d. From eq. (j5.8|) . we 
directly see that Tes diverges when Tq and Ti get closer to each other. Far from the critical region, we 
have Tcfj ^ (To — Ti)~^/^, independently of d. On the other hand, when the system enters the critical 
region, results depend on whether fluctuations are strong as is the case for 2 < d < 4 or whether one 
has mean-field criticality when d > 4. In the fluctuation-dominated regime, that is 2 < d < 4, the first 
derivative of f{p) diverges as p goes to zero. Using the known small-p expansion of / |2Hl IHll tHH], we 



in the neighbourhood of the critical point. Then we simply retain Tcs ~ (Tq — Ti) 

For 2 < d < 4, we have thus found a rather counterintuitive effect. Namely, the time scale t^s up to 
which the system evolves as if the temperature were fixed at To, increases when the final temperature 
Ti approaches the critical temperature T^. This enhancement of t^q is absent for d > 4 and must 
therefore be related to the presence of strong fluctuation effects close to Tc and below the upper critical 
dimension. 

Finally, we consider the two-time observables and the fluctuation-dissipation ratio. We have just 
seen that for r large, the kinetics of the system occurs at an infinitesimally slow rate for a long time. 
In contrast to the treatment of section 4, it is therefore useful to compare the time-evolving quantities 
with the ones when formally r = oo or, in physical terms, when the temperature remains fixed at 
Tq. Formally, the calculations can be done by analogy with section 4. First, we use the single-time 
correlator to define an effective time-evolving temperature Tes{s) from the relationship 




(5.8) 




(5.9) 



find 




limC(q,s) = 2req(To)Teff(s) 

q-^0 



(5.11) 
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We shall denote by A^°°^ the result obtained for an observable A with a time evolution governed by 
an infinitely long external time r = oo. In particular, we have \imq^QC^°°\q, s) = 2Torcq(To). As in 
section 4, it is easy to see that 



Cit,s) 



Tf 







l g{s)/g(°-\s) 



(5.12) 



where g^°°\t) is again the solution of the spherical constraint ()2.15|) but now for r = cxd and therefore 
Ti is replaced by Tq. 

We now want to find Teff(s) explicitly. Using the general expressions for the single-time correlator 
C{q,t) from section 2, the explicit form of g(t) found above and taking the limit q — 0, we obtain 



Tesis) 



Aexp{-p^s + syr^g) 
2req(ro)6 



2b 



1 + - / dti exp ( p,t, - -f ) (Ti + (To - Ti)e-*^/^) 



'eff 



Now, we consider times s which satisfy 



s T 



lA 



(To) 



(5.13) 



(5.14) 



From the explicit form of A, it is clear that this condition is the easier to satisfy the closer Tq to Ti, 
even if both are far from the critical point Tc. We then obtain 



Tesis) 



2&req(To) 
T(s) 



9h Ih p{Pc-l/r)s 

1 + -r^iPc-'e^- + -r{T, - Ti) — + O (e"^-) 



(5.15) 



where in the last step we used that r ^ Teq(To) or equivalently that Pc ^ and we recall the 
definition /(pc) = I/Tq. Therefore, we see that under the stated conditions, the effective temperature 
is equal to the temperature of the external heat bath. This nicely confirms that our definition of Tefr(s) 
is a natural one. 

Physically, this result illustrates the (expected) applicability of the adiabatic approximation to 
equal-time two-point correlators. Namely, we have seen that in the limit r ^ rcq(To), the effective 
temperature Teff(s) equals the bath temperature T{s) and the equal-time two-point correlators are 
those of a quasi-equilibrium system. 

Lastly, we examine to what extent this heuristic expectation concerning equal-time correlators carries 
over to two-time quantities. First, in an analogous manner to section 4, the two-time response function 
reads 

R^^\t,s)k^^\t,s) (5.16) 



R{t,s) = f 



9{t) 



Then, we can easily compute the fiuctuation-dissipation ratio X{t^s). Formally, it is again defined by 
()4.f)|l and we find as before 



f d 



X{t,s) \ds 



i— In C{t,s) — lnC(~)(t,s 



d 



ds 



(5.17) 



and this holds provided t — s ^ roq(To). In the just-studied adiabatic regime, we have 

C{t, s) = ZMe(^^-*^)/-o^ff C(°°)(t, s) {5Ai 



Tn 
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To r^s 2req(To)y 



1 



1 



2g req(ro) ^ req(To) (To-Ti)^^/" 
A r r Ti + (To - Ti)e-^A 



T 



(5.19) 



r To 



where the last relation is valid provided s <^ r and where the condition eq. (j5.14j) was used. 

For an equilibrium system, it is well-known that X{t, s) = 1 for all times s < t from the fluctuation- 
dissipation theorem. Here, we find an intermediate regime req(To) <C s <C r such that X(t,s) reaches 
a plateau value which is different from unity. Still, the adiabatic approximation remains valid, up to 
small corrections in s/t. However, the fact that for s ^ r the difference X{t, s) — 1 is finite (albeit it 
may numerically be small) shows that the validity of the fluctuation-dissipation relation is independent 
of, and more restrictive than, the adiabatic approximation. 

Of course, for s ^ r ^ Tcq(To) the fluctuation-dissipation theorem will be recovered. In other 
words, eq. ()5.19|) shows that the equilibration s 3> Teq(To) with respect to the internal degrees of 
freedom is already enough to reproduce the flucutation-dissipation theorem to a good approximation 
while it becomes exact when equilibration with respect to the bath, s ^ r is also achieved. 

6 Conclusions 

We have studied how the non-equilibrium kinetics of a simple ferromagnet is modified when the system 
is not instantaneously quenched to a final temperature Ti but brought into contact with a heat bath 
starting at an initial temperature Tq such that the approach towards Ti is described by a finite external 
relaxation time r, according to (jl.ip . More general temperature histories T{t) of the heat bath can 
be decomposed into sums of exponentials via Laplace transformations. Then the physically interesting 
long-time behaviour will be governed by the largest relaxation time arising in the Laplace spectrum. 
We have studied as an example the exactly solvable spherical model in order to get physical insight 
into the role of a finite value of r. We have concentrated on the high-temperature phase such that 
To > Ti > Tc, generalizing previous studies [23 1201 IZ3 123 121 with r = 0. 

Our results can be summarized as follows: 

1. If r is smaller than or at most of the same order of magnitude as the internal relaxation or 
equilibration time Tint = Teq, then for times large compared to r, Tcq the system approaches expo- 
nentially fast the known stationary state, in agreement with what was found previously for r = 



In addition, there are transient contributions to the spin-spin correlation function with spatial 
oscillations on finite length scales A. In particular, if r < Teq, we find 



2. In this case, the stationary state is indeed in thermodynamic equilibrium and is reached expo- 
nentially fast for times larger than r. We have checked this through explicit calculation of the 
fluctuation-dissipation ratio X{t, s) —>■ 1 for late times. 



[221 12SI 123 123 123- 




(6.1) 
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3. A qualitatively different behaviour may be observed in the clamped case, where r S> Teq. 

A common physical example with clamped degrees of freedom is a glass. A system in a glass 
state remains out of equilibrium as attested by the fact that the fluctuation-dissipation limit 
Xoo = lims^oo linif^oo X{t, s) ^ 1 (experimentally, this has been studied recently for the spin glass 
CdCr1.7Ino.3S4 j,19^ and a colloidal glas [20j). This means that in a glass different temperatures can 
be defined, depending on the way temperature is measured. These phenomenological properties 
characteristic of a glass can be partially reproduced in our simple model, even though we did not 
introduce disorder into the Hamiltonian. 

Specifically, we have shown, in the setting of the kinetic spherical model, that in the clamped case 

(a) the effective time-dependent temperature T^^it) as defined in eq. (|5.11|) from the equal-time 
spin-spin correlator agrees with the external bath temperature T{t) for times larger than Tgq. 

(b) the time-evolution of two-time observables is governed by the initial temperature Tq for times 
up to the new time scale (if 2 < (i < 4) 

reff-^^V^ (To-T,)-'/' (6.2) 



c 



distinct from both Tgq and r and much larger than r if Ti is close to criticality. Remarkably, 
"Tefr = i/T Tdep is the geometric mean of r and the purely internal time scale r^ep 7^ t, see 
fl5.8|l . Here the temporal evolution of the two-time observables is delayed with respect to the 
evolution of the bath temperature T{t). 

(c) there is an intermediate regime req(To) ^ s ^ r such that the system's internal degrees of 
freedom are in quasi-equilibrium while equilibrium is not yet achieved with respect to the 
temporal evolution of the bath. Then we find that the model reaches a quasi-stationary 
state, which however cannot be a quasi-equilibrium state, since the fluctuation-dissipation 
ratio 

X(M)=(l-^^^°^)" (6.3) 

takes on a value different from (albeit close to) unity in the regime req(To) ^ s ^ r. 

In distinction to glasses, the plateau value (j6.3|) is only reached in the intermediate regime 
'TeqiTo) <^ s <^ r, before the system relaxes to equilibrium for very large waiting times s ^ r. 



We stress that the clamped case arises for r large. At first sight, one might have naively expected 
that if the bath temperature T{t) is slowly changed, the system might be able to follow the evolution of 
the bath through a sequence of equilibrium states. Although our explicit results are in agreement with 
the adiabatic approximations they also show that the validity of the fluctuation-dissipation theorem 
depends on more than the validity of the adiabatic picture. This observation might be of particular 
significance for the interpretation of experimental studies on the time-evolution of cooling systems (our 
results only hold if both the initial and final temperatures are in the disordered phase). Finally, it 
remains an open question what happens when the ordered phase is entered. We hope to come back to 
this problem in the future. 

While our results were obtained in the framework of the spherical model, we expect that the quali- 
tative results, in particular on when an equilibrium state can be reached, should be of broader validity. 
Tests of this expectation in other systems would be of interest. 
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